arXiv:1503.06872v2 [cond-mat.mtrl-sci] 13Jul2015 


Spin-Orbit Torques in Two-Dimensional Rashba Ferromagnets 

A. Qaiumzadeh,^ R. A. Duine,^ and M. Titov^ 

^Radboud University, Institute for Molecules and Materials, 6525 AJ Nijmegen, The Netherlands 
^Institute for Theoretical Physics and Centre for Extreme Matter and Emergent Phenomena, 
Utrecht University, 3584 CE Utrecht, The Netherlands 
(Dated: July 14, 2015) 

Magnetization dynamics in single-domain ferromagnets can be triggered by a charge current if the 
spin-orbit coupling is sufficiently strong. We apply functional Keldysh theory to investigate spin- 
orbit torques in metallic two-dimensional Rashba ferromagnets in the presence of spin-dependent 
disorders. A reactive, anti-damping-like spin-orbit torque as well as a dissipative, field-like torque 
is calculated microscopically, to leading order in the spin-orbit interaction strength. By calculating 
the first vertex correction we show that the intrinsic anti-damping-like torque vanishes unless the 
scattering rates are spin-dependent. 
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I. INTRODUCTION 

Spin-orbitronicsii^ has attracted a lot of attention re¬ 
cently as a new subfield of spintronics^ in which the 
relativistic spin-orbit interaction (SOI) plays a central 
role. Spin-orbitronics includes generation and detec¬ 
tion of spin-polarized currents through the spin Hall 
effect the induction of non-equilibrium spin accu¬ 
mulations in non-magnetic materials through the Edel- 
stein effect the triggering of magnetization dynam¬ 
ics in single magnetic systems through spin-orbit torques 
(SOTs)f2r2i and magnonic charge pumping by means of 
inverse SOTs^ Spin-orbitronics is believed to ultimately 
enable the faster and more efficient ways of magnetiza¬ 
tion switching needed for high density data storage and 
information processing, thereby providing novel solutions 
to address the essential challenges of spintronics. In this 
paper we investigate the microscopic origin of SOTs in 
a two-dimensional (2D) metallic ferromagnet with spin- 
orbit coupling. 

The magnetization dynamics in ferromagnets is gov¬ 
erned by the seminal Landau-Lifshitz-Gilbert (LLG) 
equationji^^— 

dm dm ^ .. 

— = —jrnx Hes + aomx — + T, (1) 

where m is a unit vector along the magnetization di¬ 
rection \m\ = 1, 7 is the gyromagnetic ratio, ac is the 
Gilbert damping constant and Hgg is an effective field 
which includes the effects of the external magnetic field, 
exchange interactions, and dipole and anisotropy fields. 
The first term on the right-hand side of Eq. o describes 
the precession of the magnetization vector m around 
the effective field, while the second term describes the 
relaxation of magnetization to its equilibrium orienta¬ 
tion. Furthermore, T is a sum of different magnetization 
torques not contained in the effectieve field or damping. 

The spin-polarized current-induced magnetization dy¬ 
namics in magnetic materials arises as a result of spin 
transfer torque (STTs)ji^“— It is well known that STT 


may induce magnetization dynamics in spin-valve struc¬ 
tures, and that the exchange interaction between the 
spin-polarized current and local spins leads, e.g., to 
domain-wall motion. In uniformly magnetized single¬ 
domain systems the transfer of spin angular momen¬ 
tum from the spin-current density js to a local magne¬ 
tization is modelled by two different STT terms: i) an 
anti-damping-like (ADL) or Slonczewski in-plane torque 
T (X m X m X js, and ii) an out-of-plane field-like (FL) 
torque T (x m x js, which is typically negligible in con¬ 
ventional metallic spin valves. On the other hand, in 
ferromagnets with magnetic domains, in which spin tex¬ 
tures such as domain walls are necessarily present, the 
STT also includes reactive, T oc {js ■ V) m, and dissipa¬ 
tive, T (X m X {js ■ V)m, torquesi^^i^— 

Recently, it was demonstrated both theoretically and 
experimentally that the current-induced nonequilibrium 
spin polarization^i^ in (anti-)ferromagnets with inver¬ 
sion asymmetry may exert a so-called SOT on local¬ 
ized spins and, consequently, may lead to a non-trivial 
magnetization dynamics.— Unlike STT, the SOT phe¬ 
nomenon does not require an injection of spin current or 
the presence of spatial inhomogeneities in the magneti¬ 
zation. The magnetization switching due to SOTs may 
be achieved with current pulses as short as ^ 180 ps, 
while the critical charge current density can be as low as 
- 10^Acm-2.i2 

Quite generally Rashba SOTs can be classified as ei¬ 
ther ADL or FL torques^. The first theoretical and ex¬ 
perimental studies of SOT have demonstrated that the 
ADL SOT is proportional to the disorder strength and 
can always be regarded as a small correction to the FL 
SOT— On the other hand, in some recent experiments 
the opposite statement is made: the torques with ADL 
symmetry are more likely to be the main source of the 
observed magnetization behavior— These exper¬ 
iments are performed with ferromagnetic metals grown 
on top of a heavy metal with strong SOI and may, in 
principle, be explained by the spin Hall effect which in¬ 
duces a spin-polarized current. This spin current, in 
turn, exerts a torque on the magnetic layer via the STT 
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mechanis m^^i^^i^^ so that the ADL symmetry term plays 
the major role in the effect as discussed above. 

It is, however, a serious experimental challenge to dis¬ 
tinguish between SOT and spin-Hall STT in bilayers, 
since both torques have the same symmetryi ^^i^^ Very 
recently Kurebayashi et al^ conducted an experiment 
on the bulk of strained GaMnAs, which has an intrinsic 
crystalline asymmetry. In these experiments, the contri¬ 
bution of a possible spin-Hall-effect STT was completely 
eliminated, while sizable ADL torques were nevertheless 
detected. This provides a strong argument in favour of 
the ADL-SOT nature of the observed torque. The au¬ 
thors of Ref. attribute this torque to an intrinsic Berry 
curvature, and estimate a scattering-independent, i.e. in¬ 
trinsic, ADL-SOTi^“— This intrinsic ADL-SOT has also 
been reported by van der Bijl and Duinei^ 

In this paper we calculate both FL- and ADL-SOTs in 
a 2D Rashba ferromagnetic metal microscopically by us¬ 
ing a functional Keldysh theory approach^^ By calculat¬ 
ing the first vertex correction we show that the intrinsic 
ADL-SOT vanishes unless the impurity scattering is spin 
dependent. 

The rest of this paper is organized as follows. Section 
II introduces the model and method. In Sec. Ill we 
calculate SOTs with and without vertex corrections. We 
conclude our work in Sec. IV. 


II. MODEL HAMILTONIAN AND METHOD 


where 14 ( 1 ) strength of spin-up (down) disorder, 

and the index i labels the impurity centers R^. More 
specifically, we restrict ourselves to the gaussian limit of 
the disorder potential. 

The impurity-averaged retarded Green’s function in 
the Born approximation is given by^^— 

^k,s = + aaicTyk^ - a^ky)^ . ( 5 ) 

where = e - Sk + sM + i^s, for s =t(+) or i(-), 
<7® = (CTo-l-scrz)/2, k and e are the wavevector and energy, 
respectively, Sk = /c^/2me, and M = (A -|- Ab)/ 2. We 
have also introduced the spin-dependent scattering rate 
7 s = TTi/ninipKj^, where vq = me/27r is the density of 
states per spin for 2D electron gas, and nimp denotes the 
impurity concentration. Here we have assumed that both 
spin-orbit split bands are occupied, i.e. the Fermi energy 
is larger than magnetization splitting, ep > M. 

Following Ref. [42 we minimize the effective action on 
the Keldysh contour with respect to quantum fluctua¬ 
tions of n. This procedure gives us directly the LLG 
equation which contains torque terms in linear response 
with respect to the external field E. The effective action 
is given by A = diLp^t), where stands for the 
Keldysh contour and Lp{t) = f — H) is 

the mean-field Lagrangian. 

We further assume that we are dealing with a ferro¬ 
magnetic metal which is uniformly magnetized in the z- 
direction. Thus, we can approximate the vector n as 


We start with the 2D mean-field Hamiltonian {h = c = 

1 ), 




d?r iIjI , 


Hq + Idmp + j ■ 


( 2 ) 


where i/jl = is the Grassman coherent state 

spinor. Here, Hq is the 2D conducting ferromagnet 
Hamiltonian density in the presence of Rashba SOI)^ 


/ \ 

rir.t ^ . ( 6 ) 

In order to derive the LLG equation with torque terms 
it is sufficient to expand the effective action up to second 
order in Sn and up to first order in the vector poten¬ 
tial: = W[0(^n),A] +^,est[0('5n2),A = 0]. A 

straightforward calculation gives 


Hn = 


21 

2m, 


+ aR{(T X z) -p - ^Acr ■ nr,t - ^AbCTz, (3) 


SsOT — 


J IK^' I <J) 


where p is the 2D momentum operator, ap is the strength 
of the SOI, A and Ap are the exchange energy and 
the Zeeman splitting due to an external field in the z- 
direction, respectively, nr,t is an arbitrary unit vector 
that determines the quantization axis, and cr is the three- 
dimensional vector of Pauli matrices. 

The vector potential At = Ee~''^^ jiVL is included in 
Eq. @ to model a dc electric field in the limit —>■ 0. 
It is coupled to the current density operator, which is 
given by j = (ie/2TOe)(V — — eapcr x z, where e is 

the electron charge and me is the electron effective mass. 
Finally, the impurity potential Idmp is of the form 

= ('o W 

^ ^ i 


where Xa (a = {a;,j/}) is the response function, 



and j = -ip is the charge current density. Note that in 
the absence of SOI the term 5 'sot is 0 and only second 
order terms, S^esf^, remain. The field Sn can be split 
into the physical magnetization field dm and a quantum 
fluctuation field | as Sn^ t± = where -I- 

corresponds to the upper and — to the lower branch of 
the Keldysh contour. At first order with respect to the 
quantum component we obtain 

5sOT = 1^*1 dt'jd^r'Jd?r 7 ', 


Idmp(r) 


( 9 ) 
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FIG. 1: Feynman diagrams related to the spin-torque re¬ 
sponse function Eq. (a) undressed response function, 

and (b) the first vertex correction. The solid line corresponds 
to an electron propagator in the Born approximation, the wig- 
gly line to the coupling to vector potential and current, and 
the dashed represents a spin fluctuation. The vertical dotted 
line describes the averaging over impurity positions. 


where is the advanced component of the correla¬ 
tor, and the sum over repeated indices a is assumed. 
The LLG equation is, then, derived by minimizing the 
effective action with respect to quantum fluctuations, 
6Ses/S$, = 0. Thus, the transverse components of the 
LLG equation in the Fourier space are given by, 

+ Xo;q=0,e=0 = 0) (10) 

q=0,£ 

where represents the Fourier transformation op¬ 

erator. The functional derivative in Eq. (HOD gives 
the precession and Gilbert damping terms of the LLG 
equation,— while the second term describes the SOT. 
The dependence of Gilbert damping on SOI is second 
order in and we focus below on SOT which is of 

first order in ur. The appearance of the zero-momentum 
response function Xa;q=o,£=o in the LLG equation shows 
that the SOT is finite even for spatially uniform magne¬ 
tization, in contrast to the (non-)adiabatic STT which is 
of the first order in the gradient of magnetization. 




5Sr, 




III. CALCULATION OF SOTS 

In what follows we evaluate the spin-torque response 
function of Eq. ®, shown diagrammatically in Fig. [T1 to 
derive the SOT in the ballistic limit 7 s <C ksT, where 
ksT is the thermal energy. We calculate first the bare 
(undressed) part of the response function, depicted 
in Fig. [T^. The final result for spin torque is, then, 
obtained by adding the first vertex correction, de¬ 
picted in Fig.[T}D. Throughout the calculation we assume 
that 7 s ^ ksT ^ aakp M, where kp is the Fermi 
wavevector. The condition apkp <C M is normally ful¬ 
filled in the metallic ferromagnets of interest. Whether 
or not the condition 7 s <C kpT <C otnkp is fulfilled de¬ 
pends strongly on the sample quality. The analysis of 


spin torques in diffusive regime kpT <C 7s will require 
calculation of the full vertex correction and will be done 
elsewhere. 


A. Undressed response function 


The spin-torque response function of Eq. m without 
vertex corrections is given by 


X 


(0) 


eA 

4i 


dPk 




( 11 ) 


where = k/rrie — ctpcr x £ is the velocity vector, and 
G is the Green’s function on the Keldysh contour. From 
Eq. (fTTl) we find retarded and advanced components of 
the response function in the limit of zero frequency and 
momentum as 


4^ n^o J (27r)2 
1, 


dio dio' ■ 


fu' - U 




U -|- w' — w ± fO 

( 12 ) 


where ~ 1^® spectral func¬ 
tion and fuj = -|- 1]“^ stands for the Fermi 

distribution function. 

In the limit of weak disorder, we can decompose the 
response function into two parts: the intrinsic part Xin, 
which turns out not to depend on the scattering rate and 
describes interband transitions and the extrinsic part Xex, 
which essentially depends on disorder and corresponds to 
intraband contributions. The intrinsic part corresponds 
to the principal value integration in Eq. (HI, while 
the extrinsic part is given by the corresponding delta- 
function contribution. To leading order in ap we find 


(0)- _ 

Xin.a n »/r ^oLa, 


in,a 

eofiA 


y(0)- ^ 

Aex.a 


8M 


-Vq 


ep — M ep + M 


71 


7t 


(13a) 

(i X E)a. (13b) 


The corresponding expressions for the SOTs are the ADL 
TADl p]^ rpFL contributions, which do not take into 
account vertex corrections. 


ri(0) 

^ADL 


-^FL 


—2eaRi'om x m x {z x E), 


(14a) 


eaRApp 

M 


'ep + M 
. 7t 


ep — M 
74 . 


m X {z X Fl)(14b) 


Hence, we find that the ADL SOT in the absence of ver¬ 
tex corrections has an intrinsic origin, i.e., is disorder- 
independent. 


B. Vertex correction 

Let us now turn to the calculation of the first vertex 
correction to the spin-torque response function depicted 
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in Fig. [TJj. For the corresponding response function on 
the Keldysh contour we find 


\ f dki f dk 2 f , f , 

xGfc t,Vi^p}Gkut..t]. (15) 


(1) _ eA 

4- 


The advanced component of at zero energy and mo¬ 
mentum is, then, given by 


X 


(i)- = 


X 


eA 


-Vb 


'd'^ki 




duj duj' ■ 


fuj' - fuj 


a J (27r)^7 fl+LO — uj'—iO 

— Tr[Fl • Dki {G^^ ,^<Jb'A-k2,ujO'aG'^^^^,abAki, 

+ Gki,uj^bAk2,oj<^aAk2,oj"^bGk^^^, 


+ Aki,ujO'bGk^i^<yaGk^^^,abAki,Lj' 

+ Aki ^uj^bG k^ aAk2 ,uj' ^bG k-^ , (16) 


where the summation over the index b = {0, z} and the 
limit n —^ 0 are assumed. We have also used the nota¬ 
tions ?7o = niinp(Vt -h ^4)^/4 and r]^ = ninip(Ff - 
Using the same approximations as for the undressed part 
of the response function we obtain the intrinsic contribu¬ 
tion as 


( 1 )- ^ 


ea^A 7t + 7i p 
8M “2(7474)5 


(17) 


while the corresponding extrinsic contribution is of the 
second order in scattering rates and can be neglected. 

Thus, we obtain the FL and ADL torques in the limit 
7s "C aBkp <C M to the leading order in the SOI as 


Tfl = 


eaRVoA 

M 


ep — M ep + M 


mx (zxE^lS) 


Tadl = 


7t + 74 _ 4 


74 7t 

2eaRPQmx {mx (zx E)). (19) 


the FL- and ADL-SOTs microscopically, Eqs. m and 
(HU). The FL-SOT originates from the Fermi surface 
contribution of the response function Eq. (j8|), while 
the ADL-SOT is acquires contributions from the entire 
bands. Our main result in Eq. immediately shows 
that the intrinsic contribution to ADL-SOT is completely 
canceled in the presence of spin-independent scattering 
7t = 7i- That is, the intrinsic component of the ADL 
SOT, which originates from virtual interbranch transi¬ 
tions, is canceled by the vertex correction when weak 
spin-independent impurity scattering is taken into ac¬ 
count. Our result, therefore, explicitly elucidates the 
interplay between intrinsic and extrinsic contributions 
to ADL SOT. This result resembles the suppression of 
both spin Hall conductivity in nonmagnetic metals and 
anomalous Hall conductivity in magnetic metals, in the 
presence of spin-independent disorder.— In these ef¬ 
fects the cancelation is model dependent, and occurs for 
parabolic band dispersion and linear-in-momentum SOL 
We expect a similar scenario for intrinsic SOT. 

The existence of a Rashba effect on the interface be¬ 
tween an ultrathin ferromagnet and a heavy metal is 
the subject of intense discussion. Our results show 
that the amplitudes of the FL and ADL SOTs can be 
of the same order of magnitude depending on the rel¬ 
ative strengths of the SOI, spin-dependent scattering 
rates, and exchange interaction. Our results may qualita¬ 
tively describe Co/Pt interfaces which are characterised 
by particularly large Rashba SOI of the magnitude of 
1 eV A— Relating the strength of the Rashba cou¬ 
pling to the magnitude of the SOTs, however, would re¬ 
quire ab-initio modeling and additional experimental in¬ 
formation. Which of our results apply to more general 
models and band structures will be the subject of future 
investigation. 
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